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Contractor renormalization (CORE) is a real-space renormalization-group method to derive effective Hamil- 
tionians for microscopic models. The original CORE method is based on a real-space decomposition of the 
lattice into small blocks and the effective degrees of freedom on the lattice are tensor products of those on the 
small blocks. We present an extension of the CORE method that overcomes this restriction. Our generalization 
allows the application of CORE to derive arbitrary effective models whose Hilbert space is not just a tensor 
product of local degrees of freedom. The method is especially well suited to search for microscopic models to 
emulate low-energy exotic models and can guide the design of quantum devices. 

PACS numbers: 87.55.kd,03.67.Ac,02.70.-c,03.67.Pp 



I. INTRODUCTION 

Identifying the emergent low-energy degrees of freedom in 
a strongly-coiTelated system is a highly nontrivial problem re- 
quiring considerable physical intuition and a careful analysis 
of available experimental data |]J. The contractor renormal- 
ization (CORE) method introduced by Morningstar and We- 
instein Ij, Htl is a tool to systematically perform this task: by 
suitably selecting low-energy local degrees of freedom and 
applying a real-space renormalization procedure, one can in 
principle obtain an effective Hamiltonian which is simpler 
than the original one and therefore (ideally) more amenable 
to subsequent analytical or numerical treatment. For recent 
applications of CORE see, for example, Refs.|l|l|l|l|l|9l 

The idea behind CORE is to divide the lattice on which 
the model is defined into blocks and to retain only a small 
number of suitably chosen low-lying block eigenstates. The 
low-energy eigenstates of the full Hamiltonian defined on a 
cluster formed by two or more blocks are then projected onto 
the restricted basis formed by tensor products of the retained 
block states. By requiring that the low-energy spectrum of 
the full problem is exactly reproduced, an effective Hamilto- 
nian is obtained. The mapping onto a coarser-grained lattice 
with redefined degrees of freedom is done at the expense of 
having longer-range effective interactions. A successful ap- 
plication of the method relies on a fast decay of the effective 
interactions, which in turn depends on the correct choice of 
the effective degrees of freedom and on the particular way the 
lattice is divided into blocks, as well as on how the retained 
block eigenstates are chosen. Because physical intuition and a 
good idea of the relevant local degrees of freedom are needed 
to obtain physically sound results, we believe that this is the 
main reason why CORE has not found a more widespread use. 

The "inverse" problem of using CORE to find microscopic 
models that map well onto a desired effective low-energy 
Hamiltonian does not suffer from the aforementioned prob- 
lems: because the emergent degrees of freedom are known a 
priori and their adequacy in describing the low-energy physics 
of the device is enforced by design, the aforementioned limi- 
tations of CORE can be used to our advantage. Whenever the 



effective Hamiltonian includes sizable long-range interactions 
and/or states with a vanishing projection on the restricted basis 
are present in the low-energy spectrum, one can conclude that 
the considered microscopic model is not well approximated 
by the proposed low-energy effective model. Given current 
interest in the emulation of exotic phases via physical models 
(e.g., by using Josephson junctions or cold atomic/molecular 
gases), we expect this approach to be useful when designing 
manipulable quantum tool boxes. Finally, we note that the step 
of dividing the lattice into blocks is no longer required or even 
desirable within this context and we thus extend the method to 
models built from geometrically-constrained degrees of free- 
dom, such as quantum dimer models iflolfTTll . where the emer- 
gent degrees of freedom cannot be described in terms of tensor 
products of local states. Below we introduce an extension of 
the CORE method applicable to arbitrary basis states of the 
effective model and illustrate the application of the method 
with an array of quantum Josephson junctions lfl2ll used to 
implement a topologically protected qubit ifTsll . 



II. EXTENDED CORE METHOD 

We first review the standard CORE algorithm 
and then contrast it to the extended CORE (dubbed ENCORE) 
method proposed here. 



A. Standard CORE algorithm 

Given a Hamiltonian TL defined on a lattice C, the standard 
CORE algorithm can be described as follows: 

1. Divide the lattice C into disconnected small blocks 
B and diagonalize the Hamiltonian Ti within a single 
block, keeping M low-lying eigenstates {|^m)}i^- The 
subspace spanned by tensor products of these block 
eigenstates on a cluster C — formed by joining a num- 
ber of elementary blocks — defines the reduced Hilbert 
space within which the effective model is derived. 
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2. Diagonalize W on a cluster C consisting of N con- 
nected blocks retaining the A4 = lowest eigen- 
states {\n)}i^ with energies e„ and project them 
onto the basis formed by the tensor products of the 
block states, {|0,„j, 0,„„}}]^, forming the pro- 
jected states {|^n}}i^- 

3. Orthonormalize the obtained projected states {\ipn)}i^ 
using a Gramm-Schmidt procedure 

IV'n) = ~ X! I'^™)(V'™IV'">], (1) 

where Z„ stands for the normalization of the orthogo- 
nalized state. 

4. The range-iV renormalized Hamiltonian is then 

M 

^Y.^n\^^){^n\. (2) 
n 

By construction, this Hamiltonian has the same low- 
energy spectrum as the original one. 

5. Writing Eq. Q in terms of the tensor product states 

(/)m„)}i^' we obtain the range-A^ effective 
interactions between the blocks forming the cluster af- 
ter subtracting the previously calculated shorter-range 
interactions 

feii...ijv = - (3) 

N' = l (ji,. 

where (ii, zat/) denotes the set of all connected 
range- A^' subclusters. The effective range-A^ Hamilto- 
nian can then be written as 

= H + H + J2 h^3k + ■■■, (4) 

i {i,j,k) 

where hi is the block self-energy, hij the interaction be- 
tween nearest-neighbor blocks, hj jk a three-block cou- 
pling, etc. up to range-iV interactions. 

The successful application of the above procedure relies on a 
fast decay of the long-ranged effective interactions appearing 
in Eq. ^ and therefore one chooses the restricted set of de- 
grees of freedom by specifying the elementary blocks B and 
the retained block states {\<f)„-i)}f^ . 

B. ENCORE algorithm 

It is possible to extend the ideas presented in Sec. to 
constrained effective models — e.g., quantum dimer models, 
loop models, and string nets — for which the relevant degrees 
of freedom are no longer formed by tensor products of block 
states but, instead, by the set of configurations on a given clus- 
ter satisfying the constraints of the Hamiltonian to be emu- 
lated. We thus present an extended algorithm using alternative 
ways of selecting the restricted degrees of freedom. 



1. Choose a finite-size cluster C and build a basis 
{\<f>m)}i^ for the Hilbert space of the effective model. 
In the standard CORE method this effective basis is 
a tensor product of the relevant states on the blocks, 
whereas here it is comprised by all constrained config- 
urations on C. For example, for a quantum dimer model 
we generate all M possible dimer coverings on the clus- 
ter C. 

2. Diagonalize Tl on the cluster C, calculating the M low- 
est eigenstates with energies e„ and project 
them onto the restricted basis {\(t>m)}i^ , forming the 
projected states {\i^m)W idl- 

3. Orthonormalize by means of a Gramm-Schmidt proce- 
dure as in Eq. ([T]). 

4. The Hamiltonian within the restricted space is then 
given by Eq. (|2]i. 

5. Writing this Hamiltonian in the restricted basis 
{\4>m)}i^ we obtain the effective model 

M 

en|0m)('^m|V'n)(V'n|<^m')(<^m'|- (5) 

m,m' ,n 

It is again possible to perform a cluster expansion 
within ENCORE by using Eqs. ^ and 

Note that the above discussion is for an orthonormal restricted 
basis {|'?!>m)}i''^, such as in the example discussed in Sec. Hill 
Small changes in the procedure are required if this is not the 
case ifl^ . 



III. APPLICATION: EMULATION OF THE QUANTUM 
DIMER MODEL 

A. Array of quantum Josephson junctions 

We apply the algorithm described in Sec. lIIBI to extract the 
dimer flip amplitude t for a Josephson-junction array intro- 
duced by loffe et al. to emulate a quantum dimer model 
(QDM) lUctl on a triangular lattice. This model — first investi- 
gated by Moessner and Sondhi flTll — has the desired proper- 
ties needed to implement a topologically protected qubit and 
is given by H = Ha + + Ti.(^ with 

Ha = -tJ2[\-){//\ + \//){-\] 

(6) 

a 

with similar definitions for and Ti^ . Parallel dimers on 
the same rhombus (henceforth we refer to such configurations 
asflippable rhombi) flip with an amplitude t and interact via a 
potential strength v\ the sum runs over all rhombi with a given 
orientation. Moessner and Sondhi showed that a topologically 
ordered liquid phase exists over a finite region of the model's 
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FIG. 1: (Color online) Josephson junction array used to emulate 
the quantum dimer model on a triangular lattice (shaded lines, el- 
lipses represent the dimers; see Refs. (H and {H for details). X- 
shaped superconducting islands (thick black lines) form a kagome 
lattice with normal-state star-shaped islands (thin black lines) placed 
at the center of every hexagon of the kagome lattice. Cooper pairs 
hop between nearest-neighbor X-shaped islands with an amplitude 
given by the Josephson current Jh. A large ratio between the capac- 
itances Ci and Ch ensures an on-hexagon repulsion i?hcx to emulate 
the hard-core dimer constraint. Figure adapted from Ref. IT3I1 . 
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FIG. 2: Open-boundary clusters studied: (a) A'^ x 2 (here A'^ = 3) 
hexagon ladders; (b) ten-hexagon cluster. 



phase diagram (0.82 < v/t < 1) , somet hing confirmed in a 
number of subsequent studies llT2l[l7i[l8i 1 19ll . 

The Josephson-j unction array (JJK) can be described by the 
generalized Bose-Hubbard model 
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(7) 



where the positions of the X-shaped islands in the array are 
denoted by the indices j and k and { j, k) represents nearest- 



neighbor pairs on the kagome lattice (see Fig. [T]!. 



is the bosonic occupation number at site fj, Jh is the Joseph- 
son current between two X-shaped islands, and C is the ar- 
ray's capacitance matrix. We restrict the analysis to the case 
of hard-core bosons 113 ]. 



B. TWo-dimer flips 

In this example we focus on the off-diagonal dimer flip term 
t in Eq. (|6]l from the microscopic model [Eq. (|2]l] with the fol- 
lowing set of capacitances: C, — 1, Cx — 0.25, Ci — 2.5, 
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FIG. 3: (Color online) Nonequivalent dimer flips in the ten-hexagon 
cluster [Fig.[2jb)], comprising (a) two, (b) three, (c) four, and (d) five 
dimers. Dimer flips are represented by their associated transition 
graph: dimers (thick black lines) flip to new positions (thick light 
lines) while observing the hard-core constraint. Only the underlying 
triangular lattice of the JJK array is shown (shaded lines in Fig.[T). 
The quantity E used for gauging the validity of the mapping onto a 
QDM is defined via the multi dimer flips enclosed by dashed lines. 



and Ch = 0.5 (see Fig. [Hi. Using ENCORE these are ob- 
tained from Eq. (|5]i: the transition amplitude between dimer 
configurations | (/),„) and 10m/) (m 7^ m') is the matrix el- 
ement W*(to, to'). Our results have been obtained on the 
clusters shown in Fig.|2] 

The dominant off-diagonal term in the effective Hamilto- 
nian is the two-dimer flip with amplitude t. This process in- 
volves the creation of a virtual state with a doubly-occupied 
hexagon, with energy i^hcx, in the kagome lattice and occurs 
with amplitude t w J^/Eh^x lfl2l uM . Two-dimer flips in 
the cluster with ten hexagons are shown in Fig. Oa). Al- 
though the amplitudes for all these are ideally equal, there 
are small deviations, e.g., by the configuration of the neigh- 
boring dimers (effects of Coulomb interactions) or the open 
boundaries. All two-dimer flips depicted in Fig. [3ja) can be 
seen as being correlated and are considered individually at the 
algorithmic level. 

Figure HJa) shows results for i as a function of the Joseph- 
son coupling Jh for the various clusters and in comparison to 
second-order perturbative results. Results for the ten-hexagon 
cluster are obtained by averaging the amplitudes for the pro- 
cesses depicted in Fig.[3ja); amplitudes for the individual pro- 
cesses are shown in Fig.|4jb). The results agree up to a point 
[vertical dashed lines in Figs.|4{a) -|4{c)] where the mapping 
onto a QDM fails. This agreement is an indication that, for ca- 
pacitances and Josephson currents leading to a valid mapping, 
the low-energy physics of the JJK array is indeed described 
by a QDM with local dimer resonances. Furthermore, it also 
points to the absence of sizable finite-size effects in our re- 
sults. 
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FIG. 4: (Color online) (a) Amplitude for the two-dimer flip t in thie 
JJK array obtained from the ENCORE analysis of the finite clus- 
ters shown in Fig. |2] The (red) solid curve represents second-order 
perturbation results, (b) Results for the ten hexagon cluster are ob- 
tained as the average (triangles) of the amplitudes of the two-dimer 
processes depicted in Fig.[3ja)- (c) Added absolute values for the 
amplitudes associated to multi dimer flips (S). When S is large the 
mapping onto the QDM breaks down (vertical dashed lines). Data 
for (7. = 1, (7x = 0.25, G = 2.5, and (7h = 0.5 (adapted from 
Ref. (il). 



large amplitudes for multi dimer flips indicate that the device 
is not properly described by the effective model. We denote 
the summed absolute value of the amplitudes associated to 
these multi dimer flips by S, which are directly obtained as 
the off-diagonal matrix elements of the effective Hamiltonian 
associated to the multi dimer flips enclosed by dashed lines 
in Figs.[3lb) -[3jd). FigureUc) shows E as a function of the 
Josephson current Jh. A sudden increase in S at the same 
value of Jh for which different results for t start to deviate 
from each other in Figs. Ufa) and|4|b) indicates the breakdown 
of the mapping. The appearance of "intruder states" in the 
low-energy spectrum with negligible overlap with the hard- 
core dimer configurations also indicate the breakdown of the 
mapping. The vertical dashed line in Fig. HJc) indicates the 
point where the first intruder state appears. As Jh increases 
and charge fluctuations start to dominate, intruder states dis- 
playing multiply-occupied hexagons in the JJK array violat- 
ing the hard-core dimer constraint have their energy lowered, 
eventually causing some of the projected states {|V'm)}i^ to 
vanish. 



IV. SUMMARY 

We have presented an ENCORE algorithm suitable for con- 
strained effective models whose basis states are not simply 
tensor products of local block states. We find that CORE is 
very effective in the design of quantum devices for emulating 
exotic phases. The inadequacy of the restricted set of degrees 
of freedom in accounting for a system's low-energy behavior 
is reflected by the presence of long-range terms in the effective 
Hamiltonian obtained from CORE and is used as a criterion in 
deciding on whether successful emulation is achieved. 
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C. Multi dimer flips: Breakdown of the mapping 

Whereas a standard CORE expansion proceeds by consid- 
ering clusters comprising an increasing number of sites, an 
ENCORE expansion for the JJK array, due to the dimers' 
hard-core constraint, is performed in terms of the number of 
dimers in a cluster. The analysis of multi dimer terms can 
be used to gauge the validity of the mapping onto a QDM: 
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